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A nonintrusive laser-based measurement system has been applied for the first time in the 
HYMETS (Hypersonic Materials Environmental Test System) 400 kW arc-heated wind 
tunnel at NASA Langley Research Center. Planar laser-induced fluorescence of naturally 
occurring nitric oxide (NO) has been used to obtain instantaneous flow visualization images, 
and to make both radial and axial velocity measurements. Results are presented at selected 
facility run conditions, including some in simulated Earth atmosphere (75% nitrogen, 20% 
oxygen, 5% argon) and others in simulated Martian atmosphere (71% carbon dioxide, 24% 
nitrogen, 5% argon), for bulk enthalpies ranging from 6.5 MJ/kg to 18.4 MJ/kg. Flow 
visualization images reveal the presence of large scale unsteady flow structures, and indicate 
nitric oxide fluorescence signal over more than 70% of the core flow for bulk enthalpies 
below about 11 MJ/kg, but over less than 10% of the core flow for bulk enthalpies above 
about 16 MJ/kg. Axial velocimetry was performed using molecular tagging velocimetry 
(MTV). Axial velocities of about 3 km/s were measured along the centerline. Radial 
velocimetry was performed by scanning the wavelength of the narrowband laser and 
analyzing the resulting Doppler shift. Radial velocities of +0.5km/s were measured. 


I. Introduction 

A rcjet facilities are a vital tool in the testing and characterization of materials intended for hypersonic vehicles, 
including those designed for planetary entry. Arcjets are capable of producing flows of a higher enthalpy than 
traditional wind tunnels, and are thus better suited for evaluating and characterizing candidate materials for thermal 
protection systems. Like other types of hypersonic facilities, arcjets cannot perfectly simulate all the flow conditions 
relevant to hypersonic flight. For example, significant dissociation of freestream gases typically occurs. 
Determining how to interpret arcjet test results and then extrapolate those results to flight conditions can therefore be 
complicated. To that end, measurements of the flow conditions of an arcjet facility are needed to validate 
computational tools, to allow for reliable comparisons between simulations and test results, and ultimately to make 
possible reliable computational predictions of aerodynamic parameters and material response in flight environments 
that cannot be adequately simulated in ground test facilities. Currently, the flow parameters that can be measured 
(or calculated from measured quantities) in the Hypersonic Material Environmental Test System (HYMETS) facility 
at NASA Langley include stagnation pressure; heat flux (semi-catalytic hot-wall, fully-catalytic cold-wall, and/or 
non-catalytic cold-wall); gas mass flow rates; sonic, stagnation, and bulk specific enthalpy; and arc current, voltage, 
and power. 1 Nonintrusive measurements of additional flow parameters are therefore desired. An arcjet flow 
presents a challenging environment for making measurements as the flow is typically a high-enthalpy, low pressure, 
chemically reacting, nonequilibrium environment. Several techniques have been applied by others to arcjet flows. 
Diode laser absorption has been demonstrated for making simultaneous velocity and temperature measurements of 
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an argon arcjet plume. 2 Oxygen (O)-atom and nitrogen (N)-atom laser-induced fluorescence (LIF) have been used to 
provide temperature, velocity and species concentration (number density) measurements at a single point or along a 
line. 3,4,5 Nitric oxide (NO) planar LIF (PLIF) and O-atom LIF has previously been used in arcjets to measure 
translational temperature of O and rotational temperature of NO. 6,7 References 8 and 9 give additional descriptions 
of techniques that have been used to make nonintrusive measurements in arcjets. 

Researchers in charge of operations in HYMETS have set a goal of implementing nonintrusive diagnostics to 
obtain 1) flow visualization information, including measurements of shock standoff distance and flow uniformity 
assessments, 2) axial and radial velocity measurements in the freestream and near the test sample, 3) species 
detection and concentration measurements, including both species produced by the arcjet itself and gaseous species 
resulting from the ablation of test samples, and 4) measurements of rotational, vibrational, and electronic 
temperature. This paper presents results of the first application of NO PLIF in HYMETS and progress to varying 
degrees towards all four of these goals. Using NO PLIF, we have demonstrated that flow visualization, and radial 
and axial velocimetry can be performed over the full flowfield downstream of the nozzle exit for certain ranges of 
facility test conditions. The measurements will ultimately provide an improved understanding of the operation of 
the arc-heated facility and will also provide facility-to-facility and facility-to-flight scalings for materials tests 
occurring in the facility. 


II. Experimental Methods 


A. HYMETS Arcjet Facility 

Tests were conducted at NASA Langley Research Center in the Hypersonic Materials Environmental Test 
System (HYMETS) facility, an arcjet wind tunnel with a 400 kW power supply. When compared to other, larger 
arc-heated wind tunnel facilities, significant advantages of HYMETS include the relatively small workforce required 
to operate the facility (typically one technician), long run times (up to several hours), short down-time between runs 
and sample changes (less than one hour), and relatively low operating costs. Also, numerous optical ports offer a 
variety of views of the flow and test specimens. Figure 1 shows an overhead schematic view of the arc plasma 
generator, nozzle, and test chamber. Note that the schematic shows the retracted position of two flow probes, with 
dashed lines indicating the injected position of one of the probes. (The facility has four such probes, described in 
more detail below. When looking down the flow axis from nozzle toward the diffuser, the probes are located at the 
45°, 135°, 225°, and 315° positions with respect to the horizontal.) This schematic depicts the position of the seven 
access ports, each of which can be fitted with a 51 mm (2 inch) diameter UV-grade fused silica (quartz) window to 
provide optical access at wavelengths down to about 180 nm in the ultraviolet. Six ports are in the horizontal plane, 
one is below the centerline injected-probe position, and two are above the horizontal plane. Of the six ports in the 
horizontal plane, four are angled at approximately 45° to the flow axis and two are angled at 90° to the flow axis, 
just downstream of the leading edge of an injected probe. 

A segmented-constrictor direct-current electric arc-heater serves as an arc plasma generator. The slightly 
diverging flow issues from a convergent-divergent 8 degree half-angle Mach 5 conical copper nozzle with a 12.7 
mm (0.5 inch) diameter throat and a 63.5 mm (2.5 inch) diameter exit. Process gases consist of nitrogen (N 2 ), 
oxygen (0 2 ), carbon dioxide (C0 2 ), and argon (Ar). Test gasses are injected tangentially into the bore of the arc 
plasma generator at six discrete locations, where they are heated by a high-voltage electric-arc maintained between 
the cathode and anode to create a high temperature ionized plasma flow. The electric-arc is spin-stabilized in the arc 
plasma generator by the vortex motion of the injected test gasses. The test gasses used in the arc plasma generator 
are supplied by several compressed gas cylinders and can be custom mixed to any desired atmosphere. Adjustable 
volume percentages of N 2 and Ar are used as shield gasses near the cathode and anode, respectively, to protect the 
electrodes from rapid oxidation. 

The plasma flow from the arc plasma generator is accelerated through the nozzle and exhausted into a 0.6 m (2 
ft) diameter by 0.9 m (3 ft) long vacuum test chamber where it stagnates on one of four water-cooled 
specimen/instrumentation injection stings arranged symmetrically around the inside circumference of the test 
chamber. The flow is then captured by a collector cone with a 0.2 m (8 inch) diameter inlet plane, a 0.15 m (6 inch) 
diameter constant cross-section diffuser, and a coiled-copper tubing heat exchanger to decelerate and cool the flow. 
A two-stage, continuous-flow, high-mass-capacity, mechanical pumping system, is used to evacuate the plasma flow 
from the facility. The whole facility is cooled by a re-circulating chiller with associated booster pumps and heat 
exchangers. 

Four probes can alternately be hydraulically injected into the flow 51 mm (2 inches) downstream of the nozzle 
exit. Three of these probes typically consist of a pitot tube to measure stagnation pressure, a Gardon gauge and a 
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copper slug calorimeter to measure fully-catalytic cold-wall heat flux. The fourth probe is usually configured as 
either a Teflon® slug calorimeter to measure non-catalytic cold- wall heat flux, a silicon carbide (SiC) probe to 
measure semi-catalytic cold-wall heat flux, or a test specimen. For the results presented herein in which a probe was 
inserted into the flow, the probe used was a 25 mm diameter SiC probe. 

A more thorough description of the facility, including detailed explanations of the gas injection system, the 
instrumentation available in the facility, schematics and photographs, comparisons with other similar facilities, and 
measured free stream quantities across a wide range of flow conditions can be found in Ref 1. 



Figure 1. Schematic of the HYMETS test section. Laser sheet (shown in purple) enters test chamber through one 
of the viewing ports. A periscope (two mirrors, indicated by thick black lines) inside the test chamber then directs 
the laser sheet to the flow. Dashed lines indicate the position of one of the probes when injected into the flow. 

B. Test Conditions 

Two different gas mixtures were used for the present study. The first is used to simulate atmospheric entry 
conditions on Earth and consisted of a 75% nitrogen (N 2 ), 20% oxygen (0 2 ), 5% argon (Ar) mixture by volume. 
The second is used to simulate atmospheric entry conditions on Mars and consisted of a 71% carbon dioxide (C0 2 ), 
24% N 2 , 5% Ar mixture by volume. The total mass flow rate was varied from 76 slpm (standard liters per minute) to 
404 slpm. The arc current was varied between 100 A and 200 A. These run conditions resulted in an arc plenum 
pressure (upstream of the nozzle) of between 31 kPa and 130 kPa, and a specific bulk enthalpy between 6.5 MJ/kg 
(2,790 BTU/lbm) and 18.4 MJ/kg (7,910 BTU/lbm). (Note that the units of enthalpy are units of energy — e.g. J or 
BTU — but that the “enthalpies” referred to herein are specific enthalpies, meaning that they are actually enthalpies 
per unit mass.) Hereafter, the conditions of a given run will be referenced by the specific bulk enthalpy and by the 
test gas mixture (“Earth” or “Mars” for short). We estimate an upper bound on the average free stream static 
translational temperature to be -1,300 K (-1,900°F) for the 6.5 MJ/kg Earth condition and -1,600 K (-2,400 °F) for 
the 10.8 MJ/kg Mars condition. See section III.C.4. for an explanation of how this estimate was obtained. Table 1 
contains additional flow parameters for selected runs corresponding to cases for which specific results are shown in 
this paper. 
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Table 1. Test conditions of runs described in this paper. Asterisks indicate estimated or interpolated values 
based on measured data from similar runs. 


Run Test 

Bulk 

Arc 

Mass 

Arc 

Chamber 

Stagnation Probe 

Measurement Type 

# 

Gas Enthalpy Current 

flow 

Pressure 

Pressure 

Pressure 



MJ/kg 

Amps 

slpm 

kPa 

kPa 

kPa 


4 

Earth 

15.4 

105 

102 

34 

0.104 

0.091 

Flow Visualization 

5 

Earth 

15.4 

105 

102 

34 

0.015 

NA 

Signal Fraction 

9 

Earth 

15.9 

125 

128.6 

44 

< 0.024 

NA 

Signal Fraction 

10 

Earth 

18.4 

150 

128.6 

46 

0.127 

0.11 

Flow Visualization 

11 

Earth 

18.4 

150 

128.6 

46 

< 0.028 

NA 

Signal Fraction 

13 

Earth 

16.7 

150 

154.6 

54 

< 0.032 

NA 

Signal Fraction 

15 

Earth 

15.5 

151 

178.5 

62 

< 0.044 

NA 

Signal Fraction 

17 

Earth 

14.7 

150 

203 

69 

<0.031 

NA 

Signal Fraction 

19 

Earth 

13.7 

150 

228 

77 

< 0.029 

NA 

Signal Fraction 

21 

Earth 

12.8 

150 

254 

84 

< 0.033 

NA 

Signal Fraction 

22 

Earth 

12.8 

150 

254 

84 

0.025 

NA 

Signal Fraction 

24 

Earth 

12.0 

150 

280 

91 

0.028 

NA 

Signal Fraction 

25 

Earth 

11.4 

150 

304 

98 

0.189 

0.18 

Flow Visualization 

26 

Earth 

11.4 

150 

304 

98 

0.031 

NA 

Signal Fraction 

29 

Earth 

15.2 

200 

304 

107 

0.043 

NA 

Signal Fraction 

35 

Earth 

10.9 

174 

404 

130 

0.041 

NA 

Signal Fraction 

37 

Earth 

8.0 

125 

404 

117 

0.220 

0.21 

Flow Visualization 

39 

Earth 

8.0 

125 

404 

117 

0.039 

NA 

Signal Fraction 

44 

Earth 

6.5 

100 

403 

109 

0.039 

NA 

Signal Fraction 

77 

Mars 

14.0* 

100 

100 

38* 

0.079* 

1.90* 

Flow Visualization 

80 

Mars 

15.5 

100 

76 

31 

0.073 

1.62 

Flow Visualization 

81 

Mars 

10.8 

100 

152 

53 

0.093 

2.51 

Flow Vis and Sig Frac 

85 

Mars 

10.8 

100 

152 

53 

0.093 

2.51 

Radial Velocity 

114 

Earth 

6.5 

100 

400 

109 

0.228 

0.21 

Radial Velocity 

157 

Earth 

6.5 

100 

400 

109 

0.228 

0.21 

Axial Velocity 


C. NO PLIF Flow Visualization 

The PLIF laser system includes a tunable Nd:YAG-pumped dye laser with a Rhodamine dye mixture followed 
by a mixing crystal. Optics formed the beam into a laser sheet -50 mm wide by -0.2 mm thick (FWHM) in the 
measurement region. The laser sheet was oriented in the horizontal plane relative to the laboratory frame of 
reference and perpendicular to the axis of the primary flow. Fluorescence was imaged through the optical access 
port on the bottom of the test chamber, onto a gated, intensified CCD at a viewing angle approximately normal to 
the laser sheet. Images were acquired at 10 Hz with a lps camera gate. 

The laser was tuned to the N=13 line of the Qi branch near 225.7053 nm. (In this notation, “N” is the rotational 
quantum number of the state probed by the laser and “Q” indicates a set of transitions for which the change in 
rotational quantum number between the probed state and the laser-excited state is zero. The subscript 1 indicates that 
the parity of both the upper and lower states is positive, meaning that in both states, the electronic spin is aligned 
with the total angular momentum of the molecule.) While the fluorescence signal levels were comparable for the 
more commonly used N=l-3 lines, N=13 is a well-isolated line, which is desirable for obtaining a good Doppler- 
shift velocity measurement. 

For most arc jet conditions, we acquired 100 single shots with the sample injected. If arc jet conditions 
permitted, we also acquired some images without the sample in order to look at the core flow of the arc jet. If the 
sample is left out of the flow for too long the cooling lines in the diffuser are damaged, so obtaining images of the 
core flow with no sample was not possible at some conditions. 

Difficulties with the placement of optics (in particular, two mirrors which form a periscope inside the test 
section, shown in Fig. 1) resulted in a laser sheet which did not quite reach to the nozzle exit on the upstream edge 
of the laser sheet (although diffuse scatter off the nozzle is visible in the images, if the contrast is adjusted). This is 
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an area for improvement in later tests, and in fact was improved for the quantitative velocity measurements shown 
below. The downstream edge of the laser sheet skimmed the face of the sample. Vertically, the horizontally- 
oriented laser sheet was aligned with the center of the flow. The laser sheet was fairly uniform in intensity, although 
a few striations are noticeable in the images. 

The camera (a Princeton Instruments PIMAX-II intensified 512x512 pixel CCD camera) is effectively looking 
up through a round window port in the bottom of the test section. Since this port is directly below the sample and 
since the desired field of view is upstream of the sample, the camera is looking back toward the nozzle at an angle. 
The intensifier gate width was set to 1 microsecond with a constant gain of 250. 

D. Molecular Tagging Velocimetry (MTV) for Axial Velocity Measurements 

The NO PLIF MTV method involves writing a pattern of lines into the flowfield and observing these lines at two 
different times. The displacement of the lines is used determine the flow velocity perpendicular to the lines. To 
form a laser sheet, the collimated 226 nm beam was passed through a cylindrical lens, which focused and then 
diverged the beam, expanding it in one direction while leaving it collimated in the other. A spherical lens then 
collimated the diverging axis of the beam and focused the other axis into a thin sheet approximately 60 mm wide by 
0.5 -mm thick. To tag multiple lines of NO in the test section for a velocimetry measurement, a 50 mm long, 
LaserOptik GmbH diffusion welded lens array of 25, 1 m focal length cylindrical lenses focused the laser sheet into 
25 lines, running parallel to the model surface in the spanwise direction. The spacing of the lines was approximately 
2 mm. 

To image the tagged lines, a Cooke DiCAM-PRO camera, with an intensified 1280x1024 pixel array interline 
progressive scan CCD, was used. A 100 mm focal length, F/2.0 B. Halle Nachfl. lens was used. When used in 
double shutter mode, the camera is capable of acquiring an image pair with a minimum 500 ns delay between the 
end of the first gate and the beginning of the second. Each gate has a minimum duration of 20 ns, with delay settings 
and durations set in increments of 20 ns. A detailed discussion of the timing sequence methodology used in the NO 
PLIF experiments is provided in Ref. 10. The magnification of the images (pixels/mm) factors directly into the 
measurement of the velocity, and into the error analysis. To determine the magnification accurately, images were 
acquired of a planar surface imprinted with a regular pattern of dots that was placed in the same plane as the laser 
sheet. Use of this so-called dotcard allowed perspective distortions to be corrected as well. 10 

E. Doppler Velocimetry for Radial Velocity Measurements 

Doppler-shift-based PLIF velocimetry is an established measurement technique and has been demonstrated on 
various supersonic and hypersonic flow applications. 11 ' 15 The Doppler effect shifts the location of the spectral line 
center relative to the static vacuum center if the flow has a velocity component in the direction of the laser. This 
Doppler shift, and thus the flow velocity, can be determined from NO PLIF imaging in a variety of ways, including 
both fixed- and tunable-frequency methods. Fix-laser-frequency measurement schemes can measure velocity 
instantaneously and so are preferred when time-resolved velocity measurements are required. The process of 
scanning the laser frequency limits this method to measuring time-averaged velocity. Tuned-frequency schemes, 
however, are less susceptible to systematic error 16 and the dynamic range of the technique is not limited by the finite 
width of the spectral line or laser line as is with fixed-frequency schemes. 17 In the current study, Doppler-shift based 
velocimetry was used to obtain quantitative distributions of radial velocity for two flow conditions simulating a 10.8 
MJ/kg Mars atmosphere (Run 114) and 6.5 MJ/kg Earth atmosphere (Run 85). By scanning the laser over a small 
wavelength range, the excitation spectrum of the relatively well-isolated Qi(13) transition of NO was captured on a 
series of images using the Princeton Instruments PIMAX-2 intensified CCD camera with 512 x 512 pixel resolution. 
Specifically, the scan ranges were 3.5 pm (k L = 225.7075 — 225.7040 nm) with a scan step size of 0.05 pm and 10 
images per wavelength for the 10.8 MJ/kg Mars case (a total of 720 images acquired over 72 seconds); and 4.4 pm 
( X L = 225.7079 — 225.7038 nm) with a scan step size of 0.1 pm and 10 images per wavelength for the 6.5 MJ/kg 
Earth case (a total of 450 images acquired over 45 seconds). Post-analysis of these images revealed a shift in the 
excitation spectra, which was information used to determine absolute magnitudes of radial velocity. 

III. Analysis and Results 


A. Flow Visualization Results 

7. Single-Shot Images 

Figures 2, 3 and 4 show flow visualization results for selected enthalpy conditions for Earth and Mars, 
respectively. For each condition, three false-color single-shot images are shown. The flow is from left to right. The 
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nozzle exit is clearly visible on the left side of Fig. 2, and the face of the 1-inch diameter silicon carbide sample is 
clearly visible on the right side of the images. 



Figure 2. Location of laser sheet relative to nozzle exit and SiC probe. The intensity in this averaged image has 
been adjusted so that reflected light from the nozzle can be seen relative to the luminosity of the sample. 

The instantaneous NO PLIF images in Figs. 3 and 4 show highly non-uniform flow containing large-scale flow 
structures. These highly irregular structures appear mushroom- shaped in some images. NO fluorescence was 
observed over a large percentage of the core jet flow for lower enthalpy conditions. NO fluorescence was observed 
only on the edges of the jet flow for the highest enthalpy conditions. Background (non-LIF) luminosity was also 
observed on the face of the sample at all conditions, although the intensity of this luminosity varied significantly 
from shot to shot. 

The bow shock on the sample probe was clearly visible in many of the single-shot images, and can be seen in a 
few of the selected sample images shown in Fig. 3. The shock is evidenced by a sudden decrease in fluorescence in 
front of the sample probe. As described above, images of dotcards were acquired which allow for the determination 
of absolute scale in the images. Using an average of 100 single shots from a low-enthalpy air run where the bow 
shock in front of the stagnation probe was clearly visible (6.5 MJ/kg Earth, Run 36), the shock standoff distance was 
measured to be 9.5 + 1.1 mm. 
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Figure 4. Single-shot NO PLIF images in false color for Mars atmosphere simulations, at various enthalpies. 

2. Variation in Fluorescence Signal Intensity with Enthalpy 

One approximate indication of the concentration of NO in the flow is the magnitude of the fluorescence signal; 
while fluorescence intensity depends on many factors (including temperature, pressure, density, species mole 
fractions/quenching), in general, more fluorescence correlates with a higher concentration of NO. Figure 5 shows 
measurements of the mean signal intensity in the core region of the flow over a range of enthalpies. In order to 
calculate the mean signal intensity, averaged images were cropped to exclude all but the core of the flow. The 
boundaries defining this region of interest are depicted by a white dashed rectangle in Fig. 6. For the Earth runs, the 
SiC probe was retracted briefly to allow images to be acquired of the unperturbed free stream arcjet flow. For Mars 
conditions, the probe was not retracted, and so the images used to calculate signal fraction were of flows with the 
SiC probe injected. 

Enthalpies below about 14 MJ/kg produced detectable signal levels, whereas the images obtained at higher 
enthalpies were relatively noisy. For all conditions examined, the shot-to-shot standard deviation in the mean 
intensity level was of the same order as the mean intensity. First, the mean intensity value was calculated for each 
pixel in the measurement region. Then, the standard deviation in the intensity at each pixel in the measurement 
region was calculated. Finally, the mean value of this standard deviation was calculated for all pixels in the 
measurement region. The standard deviation relative to the mean signal may provide a qualitative indication of the 
level of unsteadiness or turbulence in the flow. Note that the signal levels were lower for Mars cases than for Earth 
cases at the same specific bulk enthalpy. One possible reason could be that the chemical pathways leading to the 
creation of NO are different for the two gas mixtures, and so the amount of NO produced with the Mars mixture is 
less than the amount produced with the Earth mixture at an equivalent enthalpy. All of the flow conditions which 
we examined in this series of tests exhibited non-uniform behavior. 
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5 10 15 20 

Specific Bulk Enthalpy (MJ/kg) 

Figure 5. Measured NO fluorescence intensity versus specific bulk enthalpy. Curve is fitted to Mean Intensity 
(Air) data only. 


3. Fraction of Core Flow with Fluorescence Signal vs. Enthalpy 

In addition to overall intensity, another metric for the usefulness of NO PLIF as a diagnostic technique in this 
facility is the percentage of the flow where fluorescence signal is observed. In order to calculate this quantity — 
hereafter called signal fraction — single shot images were cropped in the same manner as above, with the boundaries 
of the measurement region shown in Fig. 6. As in the measurement of mean signal intensity in the previous section, 
the SiC probe was retracted for Earth runs but not retracted for Mars runs. After cropping, a uniform background 
intensity value was subtracted from all the images from a given run. A threshold level of 500 counts (about 3% of 
the maximum signal intensity obtained) was then applied and the percentage of pixels with an intensity about 500 
counts was calculated. 
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33% 

13.7 MJ/kg 


15% 

14.7 MJ/kg 


12% 

15.9 MJ/kg 


9% 

16.7 MJ/kg 


7% 

18.4 MJ/kg 


Figure 7. Signal fraction measurements for Earth stimulant at various enthalpies. Percentages indicate the 
mean proportion of the core flow for each condition with nitric oxide fluorescence signal. 

Figure 7 depicts false-color and black and white cropped single-shot images for a range of enthalpies for Earth 
runs, and Fig. 8 depicts similar images for Mars runs. The false-color images show a representative single-shot for 
the given run condition (i.e. the signal fraction in the selected single-shot images is very close to the mean signal 
fraction for all images in that run). The corresponding black and white images show the effect of background 
subtraction and threshold application, with white pixels representing those with an intensity of greater than the 
threshold level of 500 counts. 


97% 
6.5 MJ/kg 


97% 
8.0 MJ/kg 


72% 
10.9 MJ/kg 


49% 
12.0 MJ/kg 


41% 
12.8 MJ/kg 
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67% 
10.8 MJ/kg 


26% 
14.0 MJ/kg 


22% 

15.5 MJ/kg 

Figure 8. Signal fraction measurements for Mars stimulant at various enthalpies. Percentages indicate the 
mean proportion of the core flow for each condition with NO LIF signal. 

Figure 9 shows a graph of the mean signal fraction as a function of enthalpy. The error bars indicate the 
standard deviation in the signal fraction for all of the single shot images used to calculate the mean signal fraction 
for each run. The exponential fitted line was fit to the Earth data; however, the Mars data appears to follow the 
same trend. Enthalpies below 10 MJ/kg had signal over almost the entire flow, but over about 12 MJ/kg, less than 
half of the flow had signal. This provides a guideline for the flow conditions where NO PLIF velocimetry 
measurements should be feasible, as well as those for which low signal fractions make single-shot molecular tagging 
velocity measurements unfeasible. This finding was in good agreement with the results of Mizuno et al. who 
reported that “total enthalpy of under 10 MJ/kg is suitable” for NO LIF in JAXA’s (the Japanese Aerospace 
Exploration Agency’s) 750kW arc heated wind tunnel. 7 One additional observation is that for two runs at identical 
conditions (the two points with an enthalpy of 12.8 MJ/kg), two different mean signal fraction values were obtained. 
The reasons for this are not entirely clear. The difference could be an indication of real variations in the arcjet flow, 
or could perhaps be a byproduct of variations in the probe laser intensity. Future measurements should include 
simultaneous measurements of laser intensity. 
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5 10 15 20 

Specific Bulk Enthalpy (MJ/kg) 

Figure 9. Fraction of core flow with NO fluorescence signal as a function of specific enthalpy. 


4. Upstream Influence of Stagnation Probe 

An unexpected finding of the flow visualization images was the effect that the presence of a stagnation probe 
had on the flow upstream of the probe. Figure 10 shows three false-color images from an 8.0 MJ/kg Earth run. 
Note that in all three images, the contrast has been enhanced over the left and right side of the images to show the 
positioning of the nozzle exit and the SiC probe (if present). The red image on the left is an average of 42 single 
shots without a probe in the flow; the green image on the right is an average of 54 single shots from the same run 
where a SiC probe is injected into the flow. The center image shows the composite of the two, highlighting the 
differences between the flow with and without the probe. With the probe in place, the NO PLIF signal extends to 
the top and bottom of the image. With the probe removed, the NO PLIF is localized closer to the core. In a purely 
supersonic flow, pressure disturbances cannot propagate upstream. However, in this wind tunnel flow, the subsonic 
boundary layer and the nearly stagnant flow outside the core of the open jet flow provide mechanisms for pressure 
disturbances to propagate upstream. Because HYMETS has a diverging conical nozzle, the flow is expected to be 
slightly diverging as it exits the nozzle. This can be seen in both the averaged image of the unperturbed core flow 
(i.e. without a probe inserted into the flow) and in the averaged image of the flow with the probe injected. 
Additionally, there appears to be relatively little NO outside the core flow in the no-probe case, but a significant 
amount of NO fluorescence can be seen in this region in the probe-in case. In the images with the probe inserted 
into the flow, the angle of divergence appears to be somewhat greater than in the no-probe case. This interpretation 
is likely misleading because the differences in the images are probably due to different amounts of NO and/or 
different quenching environments in the shear layer between the core flow and the nearly stagnant flow for the two 
different conditions. (Quenching is a process by which fluorescence intensity is reduced due to collisional energy 
transfer from excited NO molecules to other molecules. Oxygen is particularly effective at quenching NO 
fluorescence, so the additional NO present in the stagnant region in the probe-in case may displace ambient oxygen, 
thereby reducing the amount of quenching in the shear layer.) 
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Figure 10. Comparison of averaged probe-out (left) and probe-in (right) false images for Earth simulation. 

The center image is a composite of the two other images. 


B. Axial Velocity Analysis and Results 

1. Axial Velocity Uncertainty Analysis: 


The method used to process both single-shot and average velocities and associated uncertainties is similar to the 
approach outlined in Refs. 10 and 18. However, to investigate the potentially unsteady nature of the core flowfield, 
the previous method of determining the spatial uncertainty has been modified. 

In the previous estimates of single-shot spatial uncertainty, which has been documented in Ref. 10, the flow is 
assumed to remain essentially laminar in nature. This assumption led to a formulation of single-shot uncertainty 
based upon the standard deviation in the measured single-shot tagged profile shifts, which had been defined as: 


u V x ,Ax ~ 


dVx 

d{Ax) 


* u Ax — 


_ d Mk 1 


d (Ax) 


* u Ax — 


C-Atjsu) 


* tv— 1,95% * ° Ax 


( 1 ) 


In the above equations, C is the correction factor, a Ax is the standard deviation of the observed tagged profile 
shift in pixels, t^-1,95% is the student t-statistic at 95 percent confidence, A t ND is the time separation between 
sequential frames, and N being the total number of data points at a particular pixel location used to obtain o Ax . For 
this paper the calculation of correction factor, which has a dependence on the ratio of peak signal intensities between 
the initial and delayed frames, can be described by the following relation: 

1 /2 

C = a + b + c In (2) 

\^2 ,maxj \$2 ,maxj 

with the coefficients a , b , and c being 0.988, 8.993(10‘ 5 ), and -2.783(10‘ 3 ), respectively, for an effective first gate of, 
t eG 1 , of 5 ns. 

In this paper, we have obtained single-shot spatial measurement uncertainty estimates based on signal-to-noise 
considerations alone, since the signal-to-noise ratio is the largest contributor to the measurement uncertainty in 
single-shot measurements. In order to make this estimate, measurements of displacement along several tagged 
profiles were made in regions well outside of the core nozzle flowfield where the mean axial velocity was measured 
to be essentially zero for a set of 143 sequential image pairs. Additionally, over the duration of these measurements, 
the laser was tuned over the Qi(12) + Q2(20) transition from 226.025 nm to 226.020 nm. Tuning the laser away 
from the absorption peak decreased the PLIF intensity significantly, so a dependence of spatial precision uncertainty 
on fluorescence intensity could be established. By grouping together points with intensities in the first gate within 
bins, starting with points with peak intensities between 20 and 30 counts and incremented by 10-count bins up to a 
maximum of 410 counts, the standard deviation of the measured shifts corresponding to the mean peak intensity 
within each 10-count bin is obtained. Using these standard deviation values, the spatial uncertainty is calculated and 
plotted against the average intensity values for each 10-count bin, as shown in Figure 11. As the signal intensity in 
the first peak decreases, the measurement uncertainty (and the standard deviation) increases. Based upon these 
measurements, empirical curve fits of the single- shot standard deviation and uncertainty as functions of peak signal 
intensity in the first gate are obtained: 
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with a ± and ft being 1.222 and 4739.613, respectively, and a 2 and ft being 0.299 and 1238.136, respectively. 



Gate 1 Peak Signal Intensity (Counts) 

Figure 11. Determination of spatial uncertainty as a function of signal-to-noise ratio. 


To obtain an estimate of the magnitude of fluctuating axial velocity component, V x \ a relation for the 
dependence of the standard deviation of the mean velocity, must be made that separates out the signal-to-noise 
dependency identified in Fig. 11. Using the relation for standard deviation of the profile shift precision for single- 
shot measurements as a function of signal intensity put forth in Eq. (3), the standard deviation of the mean profile 
shift precision as a function of mean signal intensity, (J^x ^ max ) , is: 


^Ax(S lmax ) 


da 


Ax(S i,max ) 


dS i 


* G Si >r 


(5) 


This relation provides an upper-bound estimate of the standard deviation in the measured shift with respect to the 
mean signal at a particular pixel location. Dividing this through by the corrected time separation between the 
undelayed and delayed frames, C • At ND , and squaring the result provides the variance in the measured mean 
velocity due to random fluctuations in the mean signal intensity at the measurement location. Based upon 
observations of single-shot velocity distributions throughout each image and the corresponding signal intensities, no 
apparent correlation between the signal and measured mean velocity magnitude exists. Therefore, the variance in 
measured mean velocity due to random fluctuations in signal intensity is assumed to be independent of the variance 
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in the measured mean velocity magnitude due to random turbulent fluid mechanical fluctuations. The sum of these 
two independent variance values results in the total variance (or covariance) in the measured mean axial velocity. 

Using the signal-dependent velocity variance relation based upon Eq. (5) and the total variance in the measured 
mean axial velocity, the fluid mechanical fluctuation in the axial velocity component is: 


x j i [/— f aAx ( s i,max)\ 

V * =4 {a ^ “1 C-M N d ) 


(6) 


2. Image Pre-Processing : 

The spatial resolution for the axial and velocity experiments was measured by acquiring an image of a matrix of 
square marks separated at equal spatial intervals, known as a dotcard, mentioned above and detailed in Ref. 19. To 
correct for optical and perspective distortion of the images in these experiments, the image of the dotcard in the test 
section was acquired with the cameras and a corresponding undistorted image of the same dotcard was created with 
Adobe Acrobat software. An image registration algorithm, UnwarpJ , was then used to correct for distortion. 20 This 
software is a plug-in created for the image processing software, ImageJ , a freeware image processing program 
available from National Institutes of Health. 21 For the axial velocity measurements, once the undistorted axial 
velocity images were obtained, an undistorted background image, acquired in the absence of any fluorescence, was 
subtracted from the set. 

To improve the signal-to-noise ratio, MATLAB® was used to apply a 4-pixel radius average disk filter to the 
undelayed and delayed images. The images were then binned by 4 pixels in the span wise direction. These two steps 
improved the signal-to-noise by smoothing camera noise and consolidating the signal in regions tagged by the laser. 
However, the spatial resolution of the measurement was degraded by a factor of 4 because of this processing. 

3. Axial Velocity Processing: 

The processing of velocity data is similar to that outlined in Ref. 10. A 1 -dimensional cross-correlation method 
was used to calculate both averaged and instantaneous velocity values. In this paper, no estimations of span wise 
uncertainty have been made. However, due to the reduced spatial resolution of this experiment, the tagged profiles 
appear more closely spaced, 25 pixels peak-to-peak, than in previous analyses (Refs. 10 and 18). Additionally, the 
axial velocity magnitudes encountered in the current test are about three times larger than in prior work, resulting in 
relatively large observed shifts in the profiles between the initial and delayed frames. This presents a difficulty in 
automating the cross-correlation-based velocity measurement algorithm: a correlation window fixed about the 
tagged profile, with the same width as the peak-to-peak profile separation in the initial gate, will occasionally result 
in erroneous correlations with neighboring profiles. 

To correct for this error, a 75-pixel-wide by 1 -pixel-high window, centered about and along each tagged profile, 
was used to provide an initial estimate of profile shift. This selection of this window size ensured that 3 separate 
profiles would be encompassed, with the exception of the first and last profiles, within the initial correlation. Using 
this initial shift estimate, a reduced correlation window of 25 pixels wide was placed about the center of the initial 
profile in the undelayed frame, with an additional correlation window of the same dimension offset by the initial 
shift estimate in the delayed frame. If the peak signal intensity of the profile of initial frame was below 20 counts, 
the data was rejected.. The threshold of 20 counts was imposed to neglect regions which would otherwise return 
poor correlations and potentially skew uncertainty and fluid mechanic velocity fluctuation estimates. Signal counts 
below 20 were typically too low to produce good correlations between initial and delayed frame. 

4. Axial Velocity Results: 

Figures 12(a) and 12(b) show averaged undelayed and delayed images, respectively. These images were 
obtained by averaging the 143 pre-processed individual sequential image pairs used to compute both single-shot and 
mean axial velocities for a 400 slpm flow with a 100 amp arc current applied (6.5 MJ/kg Earth). The approximated 
time separation between the two images is 550 ns. Each line is approximately 13 pixels wide, and average 
displacements along the centerline are approximately 20 to 21 pixels. 

Figure 13(a) shows the mean axial velocities measured along the 2 nd , 1 1 th , and 20 th tagged line. The center of the 
black data points correspond to the measured mean velocities and the width of these points to the associated 
uncertainty of the mean. The portions of these profiles lying outside of the core flowfield have noticeably smaller 
(near zero) mean velocities with uncertainties, relative to the core flow region downstream of the nozzle. Along the 
11 th profile located 1.97-cm downstream of the exit plane of the nozzle, the mean velocity in this region is 
approximately 6 m/s with an uncertainty of 26 m/s. This relatively low uncertainty is attributed to the increased 
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signal-to-noise levels, persistence of signal throughout all of the sequential image pairs, and negligible fluid 
mechanic unsteadiness relative to the core nozzle flow. Within the core of the nozzle flow, both increased mean 
velocities and uncertainties are observed. For this same profile, located 1.97-cm downstream of the nozzle exit 
plane, the mean velocity and uncertainty interval within a 1-cm region centered about the axis of symmetry are 
3,073 m/s and 182 m/s, respectively. Figure 3a also provides the fluctuating velocity component, V x ', computed 
from Equation (6), and represented by the red data points. For this same 1-cm region about the centerline along the 
11 th profile, the mean value of the standard deviation, V x , is 207 m/s. This represents a fluctuating velocity that is 
approximately 6.8% of the mean velocity. 



Figure 12. Pre-processed average image from (a) undelayed and (b) delayed frame obtained during axial 
velocity measurements. The sample is seen to the right side of the image. 


Figure 13b provides a set of single-shot velocity measurement data along the same profile regions shown in Fig. 
13(a). For a 2-cm region encompassing the centerline and located 1.97-cm downstream of the nozzle exit plane, the 
measured mean single-shot velocity and uncertainty are 2,928 m/s and 760 m/s, respectively. The absence of 
continuous measurement points along the profiles displayed in Fig. 13b, and in general for the set of single-shot 
sequential images, is due to reduced signal levels within the core of the nozzle flow and intermittent presence of NO 
within this flow region. Typically within this region, only 30 to 50 percent of the total image pairs at a particular 
pixel location yield a measureable shift. Outside of the core of the nozzle flow, the approximate yield increases to 
between 90 and 100 percent. 
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Figure 13. Axial (a) average and (b) single-shot velocity profiles. The width of the black velocity points in (a) and 
(b) correspond to two times the measured uncertainty. Red data points in (a) correspond to estimated fluctuating 
velocity component. 
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Figure 14. Mean axial velocity field interpolated from MTV velocity profiles. 


Figure 14 shows the axial velocity flowfield, interpolated over the entire region encompassing the measurements 
obtained along each tagged profile shown in Fig. 12(a). The flow propagating down the centerline of the flow 
appears to have a nearly constant velocity, despite the fact that the nozzle is conical so the flow is diverging and 
continuing to expand. 

5. Radial Velocity Analysis and Results 
1. Analysis Method 

Using post processing software, Image J, a temporal distribution of the fluorescence signal was extracted from 
the image sequence at each spatial location. Since the laser scanning rate was synchronized to the camera system, a 
spectrum of the fluorescence signal versus wavelength could be obtained from the image data: each pixel in the flow 
resulting in an NO PLIF excitation spectrum. As stated in the Experimental Methods section, acquiring the data for 
a full wavelength scan at a given condition took approximately one minute (Typical spectra obtained from three 4x4 
binned regions located 4 cm downstream of the nozzle exit (PI is a region outside the core flow; P2 is a region of 
large positive radial velocity, and P3 is a region of large negative radial velocity) are shown in Fig. 15. (See Fig. 
18(a) for the location of these three regions relative to the overall flowfield.) 
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Figure 15. Three typical spectra obtained from various locations in the flow. Scatter points = experiment, solid 
line = Gaussian curve fit. Locations from which spectra are extracted are depicted as points in Fig. 18(a): PI (x = 4 
cm, y =4.4 cm), P2 (x = 4 cm, y =2.8 cm), P3 (x = 4 cm, y =-2.8 cm). 4x4 binning is used in images. 


The equation for the Gaussian curve fit shown in Fig. 15 is defined as: 


N 


p 





(7) 


where A, X L , and X D are the signal level, transition center wavelength, and the apparent transition linewidth. The 
error between the experimental data, N p>exp , and N p is defined as: 



( 8 ) 


By differentiating Eq. (8) with respect to A, X L , and X D it was possible to minimize the total error using the 
iterative Newton-Raphson method: 
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To verify that each fit variable was properly optimized, the algorithm monitored the residual error between old 
and corrected values. If the calculated error for each fit parameter was less than the specified tolerance of 1CT 10 , the 
optimization process would stop and the algorithm would proceed to the next spatial location. If the error exceeded 
the specified tolerance, the total number of iterations would double and an updated residual error would 
subsequently be calculated. For most regions of the flow, the optimization process would complete after 500 
iterations. In regions of low signal-to-noise, the algorithm would require substantially more iterations to complete 
the process. If an upper limit of 1,000,000 iterations was reached at a single location, the algorithm would discard 
that portion of the experimental data. The discarding of measurements usually occurred in regions of very low 
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signal-to-noise ratio or regions outside of the laser sheet interrogation volume (shadow regions or at the sample). By 
filtering out regions with low signal-to-noise (i.e. A < 25) prior to optimization, faster processing was achieved. 

After obtaining the values of X L and X D it is possible to calculate quantitative values of radial velocity. Since 
the laser sheet was projected normal to the axial hypersonic flow produced by the facility, any radial component of 
velocity will cause an apparent shift in the transition center frequency, v C9 relative to the laser center frequency, v L . 
This shift is termed the Doppler shift and is defined as: 


Av n =v-v, 


( 12 ) 


The wavelength, A, and Doppler shift in wavelength, A, are related to wavenumber, v, and Doppler shift in 
wavenumber, v, through the following equations: 
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The Doppler shift is related to the component of radial velocity, V n and the speed of light, c, by: 

r v, 


(15) 


The measured shift in the excitation spectra were assumed to be due entirely to the Doppler shift, thereby 
ignoring the collisional shift (also known as pressure shift), which is significant at higher pressures. Errors resulting 
from this assumption are considered below. Also, in the analysis, symmetry was invoked only to determine where 
the centerline frequency reference was for stationary molecules. By symmetry we assumed that the radial flow 
velocity in the outermost regions of the flow were equal in magnitude but opposite in sign. However, the Doppler 
shift was observed to be equal in these two regions. Thus the velocity in both of these regions must be zero, 
providing a convenient reference for zero Doppler shift. 


2. Uncertainty Analysis 

The main sources of random error are due to randomness in the data and low signal-to-noise ratio. As the 

tolerance in the optimization algorithm was set to 1CT 10 , the corresponding uncertainty in the velocity is 
approximately ±10‘ 5 m/s, which is negligible. The uncertainty due to low signal-to-noise is demonstrated by looking 
at the effect of binning (Fig. 16). The effect of binning is to increase the total signal-to-noise ratio, but also to 

decrease the spatial resolution of the velocity measurement. 



Figure 16. Effect of binning on the radial distribution of radial velocity (10.8 MJ/kg Mars, Run 114). 
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It was determined that a 4X4 binning provided adequate spatial resolution to capture any small scale structures 
independent of the measurement technique, but also provided sufficient integration of the signal to reduce random 
noise observed in the unbinned data of Fig. 16. Regions in the flow at the centerline and near the top and bottom 
wall surfaces have a near zero radial velocity component. By analyzing the fluctuation in velocity in those regions, 
an estimate for the uncertainty due to random noise and flow fluctuations could be determined. Near the top and 
bottom wall surfaces the uncertainty was ±16.8 m/s (95% confidence). A more conservative estimate of the error 
due to random noise and flow fluctuations could be determined from analyzing the variation in radial velocity near 
the flow centerline. This uncertainty was ±26 m/s. 

Sources of systematic error in the measurement include laser beam attenuation, wavelength linearity of laser, and 
collisional shift. Due to Beer’s law of absorption, the laser intensity decreases approximately exponentially as the 
laser passes through the flow. The black curve in Fig. 17a shows the measured radial distribution of signal level (A) 
at a streamwise position of 4 cm downstream of the nozzle exit. By assuming that the absorption coefficient is 
constant along the path of the light and by imposing symmetry on the fluorescence intensity profile (in the absence 
of laser absorption), it is possible to calculate a radial profile of laser intensity (green curve in Fig. 17a). Correcting 
the black curve by the attenuated laser beam recovers a nearly symmetric red fluorescence profile (red curve in Fig 
17a). The attenuation in laser intensity will be more significant near the center of the absorption transition. Doppler 
shifted spectra near the bottom of the image will therefore have a different laser intensity that varies depending on 
the Doppler shift (and resulting absorption) of the gas above. This error is simulated in Fig 17b. After applying a 
simulated laser intensity profile to the experimental data and re-fitting the shifted spectra (shown as dashed lines in 
Fig 17b), a systematic error in A, c can be calculated (see red and blue curves in Fig. 17b). Laser beam attenuation acts 
to artificially amplify v D in regions of high absorption and large radial velocity, which results in an over prediction 
in the magnitude of radial velocity. The maximum uncertainty due to laser beam attenuation was measured to be ±56 
m/s. 


(a) Laser Beam Attenuation 



(b) 



Figure 17. Effect of laser beam attenuation on the measurement of Doppler shift. 


There is also an uncertainty in the radial velocity measurement due to non-linearity in the laser scanning drive. 
Such an error would be caused by the laser’ s drive software not indicating the actual change in wavelength (absolute 
wavelength inaccuracy does not lead to an error in the current measurement). Based on the scan linearity 
measurements obtained from the manufacturer, a worst case error of approximately ±0.0045 cm' 1 in v D per wave 
number scanned was measured. This is equivalent to a systematic error in radial velocity of approximately 
±0.009(V r ). Thus, the maximum error in velocity due to non-linearity of the laser scanning is ±4.5 m/s. 

The collisional shift also contributes to a systematic error in radial velocity measurements. Though this is a 
chemically reacting flow for which a perfect gas analysis does not strictly apply, such an analysis can be useful for 
estimating static conditions in the jet for the purposes of this uncertainty analysis. For an isentropically expanded 
Mach 5 jet having a stagnation pressure (arc pressure) of 1.1 atm, the pressure shift at the nozzle exit is 
approximately 0.0006 cm-1, which corresponds to systematic error of just 1 m/s. This error was partially mitigated 
by the method of which the center wavelength, X L , was calculated. Since A, L is calculated from regions in the flow 
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without a radial velocity component, any pressure shift resulting from a large uniform pressure field would be 
accounted for. 

In summary, the total uncertainty in the flow depends somewhat on the spatial location. The maximum 
uncertainty measured in the flow is ± 62 m/s, which includes both random and systematic errors. This error is 
dominated by the laser beam attenuation. Table 2 shows the contribution of errors at three representative locations 
in the flow. In future experiments, weaker NO transitions could be probed to minimize the error caused by 
absorption. 


Table 2 Contribution of random and systematic errors at three positions (P1-P3 located at x = 4 cm) 


Uncertainty Type 

PI (y = 4.4 cm) 

P2 (y = 2.8 cm) 

P3 (y = -2.8 cm) 

Random Noise (m/s) 

25.9 

25.9 

25.9 

Laser Beam Attenuation (m/s) 

0.1 

14.6 

55.9 

Laser Non-Linearity (m/s) 

0.0 

3.0 

4.5 

Collisional Shift (m/s) 

1.1 

1.1 

1.1 

Total Uncertainty (m/s) 

25.9 

29.9 

61.8 


3. Radial Velocity Results 

Fig. 18(a) shows a vector plot for the 6.5 MJ/kg Earth condition (Run 114) overlaid with contours of radial 
velocity. The streamwise component of velocity was determined from the MTV technique (Run 157). Fig. 18(b) 
shows radial distributions of radial velocity at four different streamwise locations labeled L1-L4. The position of 
lines L1-L4 is shown as white dashed lines in Fig. 18(a). The white points, labeled P1-P3 are the positions where the 
spectra shown in Fig. 15 were obtained. Due to a varying systematic error, the magnitude of the error bars in Fig. 
18(b) also vary with radial and streamwise position. Due to the systematic nature of these errors, the error bars are 
also positioned asymmetrically to reflect the direction that laser beam attenuation and laser wavelength non-linearity 
affect radial velocity. The regions of Fig. 18(a) not containing measurements correspond to regions of the flow 
where the optimization algorithm discarded data due to low signal. 


V ref = 3000 m/s 




Streamwise Position (cm) 



Figure 18. Map (a) and profile (b) distributions of radial velocity for 6.5 MJ/kg Earth condition (Run 114). 

White dashed lines shown on contour map are the locations for the quantitative distributions. White points on 
contour map are locations for sample spectra fits shown in Fig. 14. 
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Figure 19 shows a comparison between the velocity field measured from the 6.5 MJ/kg Earth condition (black 
vectors) to a theoretical velocity field (red vectors) assuming a perfectly expanded jet with a nozzle half-angle of 8 
degrees. To construct the field, a point source was located upstream of the nozzle, the location of which was 
determined by projecting the nozzle walls back to a single point. In the theoretical velocity field, the total velocity 
magnitude was considered constant downstream of the nozzle at a value of 2994 m/s, which corresponds to the 
average centerline streamwise velocity in the experimental data. Although qualitative, there is close agreement 
between the two cases, which confirms the point- source nature of the flow produced by the facility and gives overall 
confidence to both the Doppler-shift based velocimetry and MTV techniques. 



Figure 19. Comparison of experimental velocity field (black vectors) to an ideal velocity field assuming a 
perfectly diverging field of constant velocity magnitude with a unique point source (white vectors). 

Experimental data correspond to the 6.5 MJ/kg Earth condition (Run 114 and Run 157). 

Figure 20 shows a comparison of maps of measured radial velocity for facility runs corresponding to the 6.5 
MJ/kg Earth condition (Run 114) and the 10.8 MJ/kg Mars condition (Run 85). Although the overall magnitudes in 
radial velocity are similar in the top portion of the image, the magnitude of radial velocity measured for the 6.5 
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MJ/kg Earth condition is 25% larger in the lower portion of the image. For the 10.8 MJ/kg Mars condition, the 
magnitudes in radial velocity are similar for both the top and bottom portions of the image. This may indicate that 


the systematic error in radial velocity due to laser beam attenuation is larger for the 6.5 MJ/kg Earth condition 
compared to the 10.8 MJ/kg Mars condition. 
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Figure 20. Comparison of radial distribution of velocity shown as a contour plot for the 6.5 MJ/kg Earth 
condition (Run 114) and 10.8 MJ/kg Earth condition (Run 85). 


4. Free Stream Static Temperature Estimation 

From measurements of the spectral width, Av D , of the Doppler-broadened Gaussian distribution such as shown in 
Fig. 15, it is possible to extract the translational temperature, T, through the following relationship: 


1 8kQn(2))N A 
\j (MW N0 )c 


(16) 

where k, N A , MW N0 are the Boltzmann constant, Avogadro’s number, and the molecular weight of NO, respectively. 
The average translational temperature in the core of the jet at the nozzle exit was measured to be -1,300 K (-1,900 
°F) for the 6.5 MJ/kg Earth condition and 1,600 K (2,400 °F) for the 10.8 MJ/kg Mars condition. When determining 
Av d in the fitting process, collisional broadening was neglected. Additionally, we assumed that the laser’s line width 
of 0.07 cm' 1 added in quadrature with Av D to fit the measured linewidth. The effect of laser beam absorption and 
saturation are systematic errors that both artificially broaden the transition. Therefore, the actual translational 
temperature is likely cooler than measured here if saturation and absorption are significant. Thus, the above reported 
temperatuer is an estimation of the upper limit of the average free stream static translational temperature. 


IV. Conclusion 

NO PEIF has been successfully applied to the HYMETS facility at NASA Langley Research Center for the first 
time. Flow visualization and velocity measurements have helped to characterize previously uncertain flow 
parameters by providing both qualitative and quantitative temporally and spatially resolved information about the 
arcjet flow. We believe that these are the first NO PLIF based axial and radial velocity measurements in an arc- 
heated facility. The techniques demonstrated in this paper are expected to be applicable in arc -jet flows having 
enthalpies of less than 10 MJ/kg, where most images exhibit strong fluorescence. At higher enthalpies, O-atom and 
N-atom LIF can be used to determine flow properties, although these are point measurements as opposed to planar 
measurements. 
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